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Abstract 

We construct an exactly solvable example of Sturmian bound states which exist in 
the absence of any confining potential. Their origin is topological - these states are 
found to live on certain "knotted" contours C^ N ' of complexified coordinates. 
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1 Introduction 



A physical framework and motivation of our forthcoming considerations lies partially 
(though not only) in the standard radial Schrodinger equation 

- ^ m + m + a m = e m (i) 

where the radial coordinate £ runs over the half-axis R + and where the standard 
Dirichlet boundary conditions are usually imposed at £ = oo and in the origin 
(at £ = 0, with the well known exceptions for strongly singular V(l; 2 ) [1]). In 
phenomenological setting, such an ordinary differential equation is usually obtained 
from a "realistic", spherically symmetric D— dimensional single-particle Hamiltonian 
H = — A + A y(|a;| 2 ) acting in the most common representation L 2 (R D ) of the Hilbert 
space of bound states. With £ = (D — 3)/2 + m in the m— th partial wave, one has 
to distinguish between D = 1 and D > 1 [2J. At D = 1 the situation is exceptional 
and a due care is needed when one tries to work, purely formally, with m = (for 
the even-parity states) and m = 1 (for the odd-parity states). At all the higher 
dimensions D > 2, the correspondence between H and eq. ([1]) is more standard and 

the sequence of the angular-momentum indices becomes infinite, m — 0, 1, In 

practice, one usually works with a fixed strength A = 1 of the interaction and studies 
the spectrum of the bound-state energies E n , n = 0, 1, . . .. Alternatively, one can 
choose and fix the energy (say, E = 1) and compute the related eigencouplings A n , 
n = 0, 1, . . . which correspond to the normalizable solutions VKO called Sturmians 
(this will also be our choice in what follows). 

An immediate mathematical inspiration of our present note can be traced back 
to the 1993 paper by Buslaev and Grecchi [2] who complexified, purely formally, the 
variable £ in eq. (pQ) (for more details cf. section [2] below). From the historical point 
of view it proved a bit unfortunate that the Buslaev's and Grecchi's considerations 
did not attract too much attention. It took further five years before Bender, Milton 
and Boettcher [5] returned to the subject, revealed and emphasized its formal appeal 
and persuaded many physicists about many phenomenological potentialities hidden 
in similar models. In this setting, our present brief note can be read as a part and 
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a continuation of the deeper analysis of the redefined models ([T]) using complex £ 
which, strictly speaking, shouldn't be called a "particle coordinate" anymore [I]. 

We shall study the maximally simplified, analytically solvable version of the 
Schrodinger differential equation without any interaction, 

-^m+^^-m=mo- (2) 

For compensation, the complexification of coordinates will be assumed more sophis- 
ticated than usual. In the formal definition of £ G C^ N \ the complex contours C^ N ' 
will be specified as highly unusual and topologically nontrivial (cf. section [3] be- 
low). As a consequence, we shall be able to obtain bound states by imposing the 
corresponding more or less standard complexified asymptotic boundary conditions 
(cf. a broader context outlined in refs. [3] and/or [5j (6J [7]). The simplicity of the 
dynamics encoded in eq. (J2J) will enable us to construct our bound-state solutions 
in closed analytic form (cf. section Hj). A more detailed discussion concerning the 
interpretation and perspectives of applicability of our G L 2 (C^) will be added 
in section [5] and in a brief summary. 

Marginally, let us note that sometimes, one could need a slight extension of the 
scope of our model (jSJ) beyond its purely kinematical version. This can be easily 
achieved by an addition of a trivial potential V(r) = ^/r 2 and by the subsequent 
redefinition of the effective i in (j2J), 

£(£+l)= 7+ ( m + :^_^) (m + E^j , m = 0,l,.... (3) 

In this way one can treat I = £(7) as a continuous, not necessarily just a (half)integer 
real parameter. 

2 The Buslaev's and Grecchi's model as a guide 
2.1 Isospectral Hamiltonians 

Let us briefly return to the Buslaev's and Grecchi's paper [2] where a constant 
shift e > has been used to define the following straight line of "unmeasurable" 
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complexified coordinates, 

COBG) = {£ = x-ie\e>0,x<ElR}. (4) 

A very specific anharmonic-oscillator potential has further been chosen as acting 
along C^ BG \ Under the most common Dirichlet asymptotic boundary conditions one 
reveals that with -0[£(±oo)] = we have (jBG) (x) = t/>[£(x)] G L 2 (C {bg) ) obtainable 
from the differential equation 

0< BG >(a;) = O. (5) 

This is a non-Schrodinger, non-self adjoint eigenvalue problem with VT— symmetry 
defined in terms of the spatial reflection V and temporal reflection T and exhibited 
by the Hamiltonian H^ BG) [2]. 

After a "naive" choice of the Hilbert space n {ori9inal) = L 2 (C {BG) ) the Buslaev's 
and Grecchi's Hamiltonian H^ BG ^ proves manifestly non-Hermitian (and, hence, ap- 
parently "unphysical" ) . Fortunately, one of the main results of ref. [2] tells us that 
jj(bg) proves isospectral to another operator 

h (BG) = nH (BG) n -l (6) 

which happens to be self-adjoint and, hence, physical. This observation settled the 
questions of physics beyond BG model and re-established the correct probabilistic 
interpretation of all the observables in the system in question. 

Several papers (cf., e.g., [8] or [9]) re-analysed the Buslaev's and Grecchi's conclu- 
sions recently. This partially motivated also our forthcoming considerations. One of 
our reasons was that for the model H^ BG ^ it was trivial to guarantee, by construction, 
that the Hamiltonian h^ BG ^ becomes self-adjoint in its own Hilbert space i-[(p h y stcal ) . 
The challenge of a search for some other simple models was imminent. 

2.2 General formalism and an amended Dirac's notation 

On a formal level needed in our forthcoming considerations one should refer to the 
review paper [10] where the authors emphasized that h^ BG ^ in ([6]) can be self-adjoint 
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dx 2 [£(x)] 2 



±^l + XV^{[ax)} 2 }-E^ 



(in the Dirac's transposition-plus-complex-conjugation sense, i.e., 

h (BG) = ^(BGjy) 

only if ij( BG ) is quasi- Hermitian (i.e., only if (tf( BG )) = O 1 where, in our 

present notation, = f^fi). From such a point of view, the Buslaev's and Grecchi's 
original choice of their very specific anharmonic-oscillator model can be interpreted 
as a "mixed blessing". On the negative side, the narrow-minded results of ref. [2] 
did not prove too inspiring. In fact, they looked so exceptional that the physics 
community accepted them as a mere mathematical curiosity. On the positive side, 
the tractability of the model seems to have opened new perspectives. 

The point is that in principle, the lower-case "correct" Hamiltonian operator 
can be interpreted as acting in another, different Hilbert space 

^physical) _ Th 

the Hamiltonian h^ BG ^ can be, in general, very different from its original upper-case 
representation introduced as acting in a "tentative" , unitarily non-equivalent Hilbert 
space 7^(°"9 ma . Exceedingly complicated versions of the "physical" h^ BG ^ may be 
encountered in some realistic models, e.g., in nuclear physics [TO] . 

In the standard Dirac's notation all the elements of the original vector space 
and of its dual may be treated and denoted as the usual kets |^) and 

bras (^1, respectively. After the change of the spaces 7-[(° rt 9 mal ) — > ^{physical) ^ - IS 
necessary to keep the trace of the changes in order to avoid the possible ambiguity of 
the notation. Thus, the elements of the physical Hilbert space will be denoted here 
by the specific, curly ket symbols G ^{(phystcai) w j 1 j] e j n fa e cl ua i space of linear 
functionals [EE] we shall write G (i^hysicai)^ 

Puzzling as it may seem at the first sight, our emphasis on the difference between 
the spaces of kets |tf y G H {vhysical) and e H {ori9inal) did in fact play a key role in 
some misunderstandings which appeared in the current literature [12J. Paradoxically, 
the most natural and transparent resolution of the whole puzzle is virtually trivial. 
According to our recent proposal [T3] one can simply add an auxiliary, third Hilbert 
space fi( third ) exhibiting the following properties: 

• as a vector space without inner product, the set 7-[( third ) coincides with / }{( on 3 mal ) ; 
i.e., we may write |*) e H (physical) as well; 
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the spaces of duals (or, if you wish, linear functionals) are different, i.e., 

^j^(original) \ ' ^(original) '^(original) ^J^(third)\ ■y (third) y^ithird) . 

for the auxiliary, innovated conjugation * and functionals ((^| G (^H^ third ^ we 
have to postulate the defining relation 



«¥|:= (*|0 = mf * (l*}) f 



e = fit n . 



(7) 



As long as the simultaneous use of both the conjugations would almost certainly 
lead to dangerous confusions, we shall always employ just the Dirac's transposition- 
plus-complex-conjugation one here. This means that 7~(°"9 ma or 'X{p h v stcal ) w iH be 
both characterized by the same single-cross superscripts ' and by the usual bra-ket 
correspondence. In contrast, the double cross * will not be used at all. Thus, we 
shall always treat the space h^- thlrd )^j and its double-bra elements ((^/\ ^ (^\ as 
mere abbreviations. 

All these conventions are summarized in Table [TJ They immediately imply that 
-< ipW {{^W) so that the spaces f{(p h y sical ) anc l j-{(third) are ^ construction, 
unitarily equivalent. Any one of them may be employed as physical, therefore. Of 
course, the same language and physical interpretation is applicable not only to the 
BG model but also to all the models with real spectra (including the quantum knots 
to be described below) which may only look non-Hermitian due to the naive initial 
choice of the "wrong" inner product in U iori9inal) . 



Table 1: The triplet of Hilbert spaces in quantum mechanics 



Hilbert space 


element 


dual 


inner product 


Hamiltonian 


'^(original) 




m = m j 


<v#'> 


H (BG) ^ ^(BG)^ 


'^(physical) 




-<ti)\ = {ifj\tf 


-<ii)\ii)'>- 


h (BG) = ^(BG)jt 


^(third) 






«V#'> 


H (BG) _ f H (BG)\* 



6 




Figure 1: Sample of the curve with N — 1. 

3 Integration paths C^ N > 

Let us now return to our ordinary linear differential eq. (T5]) of second order, the 
general solution of which can always be expressed as a superposition of some of its 
two linearly independent components. In particular, in a small complex vicinity of 
the origin we may write 

^(0 = c + ^ (+) (0 + C-^ ( - ) (0 (8) 

where 

r ^ + \C) = C, e+1 + corrections , ip^~\r) = + corrections , |£| <S 1 . (9) 

In the asymptotic domain we shall prefer another option with k = \fE in 

m = c^{i)+c 2 ^\o (io) 

where 

^ (0 = ex P (=ti K i + corrections) + corrections , |£| ^> 1 . (11) 
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Figure 2: Sample of the curve with N = 2. 

In between these two extremes, our differential eq. ([2]) is smooth and analytic so that 
we may expect that all its solutions are locally analytic. 

In the vicinity of the origin £ = our centrifugal pole with its real parameter 
£ dominates our eq. (|2J). Once i is assumed irrational, both the components of our 
wave functions (as well as their arbitrary superpositions) would behave, globally, as 
multivalued analytic functions defined on a certain multisheeted Riemann surface TZ. 
In the other words, our wave functions would possess a logarithmic branch point in 
the origin, i.e. a branch point with an infinite number of Riemann sheets connected 
at this point [T3] . 

Separately, one should study the simplified models with the rational is which 
correspond to the presence of an algebraic branch point at £ = 0. A finite number 
of sheets jH] would be connected there. In the simplest possible scenario of such a 
type we may take £(£ + 1) = with either £ = — 1 or £ = 0. In such a setting, eq. ([8]) 
just separates ^>(£) into its even and odd parts so that the Riemann surface itself 
remains trivial, 1Z = (C. 

In the generic case of a multisheeted 1Z we intend to show that the asymptotically 
free form of our differential eq. d2J) with the independent solutions (ITT]) can generate 
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bound states. One must exclude, of course, the contours running, asymptotically, 
along the real line of £ since, in such a case, both our independent solutions ip^ 1 ' 2 ^ (£) 
remain oscillatory and non-localizable. The same exclusion applies to the parallel, 
horizontal lines C^ BG ^ in the complex plane of £. In the search for bound states, both 
the "initial" and "final" asymptotic branches of our integration paths C^ N ' must have 
the specific straight-line form £ = ±|s| e ltp with a non-integer ratio p/ir. Thus, we 
may divide the asymptotic part of the complex Riemann surface of £ G TZ into the 
sequence of asymptotic sectors 

<S = {£ = -i Q<* v | Q > 1 , if e (-tt/2, vr/2)} . (12) 

S±k = {£ = -ie ±ifc >e i ^|f ? > 1, <p G (-7r/2,7r/2)}, A; = 1,2,.... (13) 

We are now prepared to define the integration contours C^ N \ For the sake of con- 
venience we shall set all their "left" asymptotic branches in the same sector 
S and specify £ = (s + s ) (1 + is) where s G (— oo, — s ), e > and s > 0. The 
subsequent middle part of C^ N ' must make counterclockwise rotations around the 
origin inside 1Z while s G (— sq,so). Finally, the "outcoming" or "right" asymptotic 
branch of our integration contour C^ N ' with s G (so, oo) must lie in another sector 
S2N of 7?., i.e., in the Riemann sheet where the requirement of VT— symmetry [6] 
forces us to set £ = (s — s ) (1 — ie). 

4 Bound states along nontrivial paths 

From now on we shall assume that the integration contour is fixed and that the 
variability of £ is confined to it. In this spirit we also adapt our notation writing 
£ = £( s ) = ^ G C^ N \ Our illustrative Figures [T] - [3] sample the choice of N — 1, 
N = 2 and iV = 3, respectively. 

It remains for us to impose the asymptotic boundary conditions requiring that 
our wave functions vanish at s — > ±00. As long as our integration path C^ N ' per- 
forms counterclockwise rotations around the origin, this form of the asymptotic 
boundary conditions will already guarantee the normalizability of our bound-state 
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Re r 



Figure 3: Sample of the curve with N — 3. 



wave functions ip(r) G L 2 (C^J (cf. the similar situation encountered in the models 
with confining potentials [2j [5j [15]). 

In our bound-state problem (j2J) with the (by assumption, real) E = k 2 we may 
set z = Kr and V( r ) — a/^V 9 ^)- This reduces eq. (j2J) to the Bessel differential 
equation with the pair of the two well known independent special-function (say, 
Hankel- function [16]) solutions which may be inserted in our ansatz 



^(r) = c 1 V^Hi 1) (Kr) + c 2 V^Hl 2 \Kr), v = £+1/2. 



(14) 



At |argz| < 7r and Rez/ > —1/2, the asymptotics of its components are given by the 
respective formulae 8.451.3 and 8.451.4 of ref. [T6j . 

/ vr(2z/ + l)^ 



V /ffli 1) (z)=ex ] 
yf^ 2 )(,)=exp [-i( 



7r(2z/ + r 



1 + 



2iz 

2 1/4 



2iz 



This implies that inside the even-subscripted sectors S 2 k our ansatz f[T4|) combines 
the asymptotically growing (and, hence, unphysical) component Hjp(z) with the 
asymptotically vanishing and normalizable, physical component H^'(z). Vice versa, 
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in all the odd-subscripted sectors <S>2fc+i we would have to eliminate, in principle, the 
asymptotically growing H^\z) and to keep the asymptotically vanishing Hjp(z). 

We may start our discussion of the existence of the localized bound states from 
the straight-line contour C = C^ BG > = which is all contained in the zeroth sector 
Sq. This immediately implies that with Imr <C —1, the asymptotically vanishing 
solution ip^(r) = H^\kt) remains unconstrained at all the real k. Obviously, 
the spectrum remains non-empty and bounded from below. This means that the 
low-lying states remain stable with respect to a random perturbation. A less usual 
feature of such a model is that its energies densely cover all the real half-line M + . This 
feature is fairly interesting per se, although a more detailed analysis of its possible 
physical consequences lies already beyond the scope of our present brief note. 

Our eigenvalue problem becomes not too much more complicated when we turn 
attention to the spiral- or knot-shaped integration contours with N > 0. In such 
a case, fortunately, the exact solvability of our differential equation enables us to re- 
write ansatz (ITU]) in its fully explicit form which remains analytic on all our Riemann 
surface 71. Once we choose our "left" asymptotic sector as So, the "left" physical 
boundary condition fixes and determines the acceptable solution on the initial sheet, 

^(r) = Cy /fHl 2) (Kr), rE S . (15) 

After the N counterclockwise turns of our integration path & N > around the origin 
this solution gets transformed in accordance with formula 8.476.7 of ref. [16] which 
plays a key role also in some other solvable models |17j . 

/#> (ze™) = Sin(1+m)7n/ H?\z) + H?\z) . (16) 

Here we have to set m = 2N. This means that the existence of a bound state will 
be guaranteed whenever we satisfy the "right" physical boundary condition, i.e., 
whenever we satisfy the elementary requirement of the absence of the unphysical 
component HjP(z) in the right-hand side of eq. (fT6|) . 

The latter requirement is equivalent to the doublet of conditions 

2Nu = integer , v ^ integer . (17) 
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This means that at any fixed and positive value of the energy E = k? and at any 
fixed winding number N = 1, 2, . . ., our present quantum- knot model generates the 
series of the bound states at certain irregular sequence of angular momenta avoiding 
some "forbidden" values, 

£= ^27\T~' U 1 U / 2X. \X.(i.X (18) 

These bound states exist and have the analytically continued Hankel-function form 
ffT5j) if and only if the kinematical input represented by the angular momenta i is 
restricted to the subset represented by formula (fl8|) . 

Our construction is completed. Once we restrict our attention to the purely 
kinematic model with 7 = 0, we can summarize that at the odd dimensions D = 2p+l 
giving £ = n+p — 3/2 we may choose any index n and verify that formula ffl8l) can be 
read as a definition of the integer quantity M = (2n+2p— 1) N which is not forbidden. 
At the even dimensions D = 2p we equally easily verify that the resulting M is always 
forbidden so that our quantum-knot bound states do not exist at V(r) = at all. 

The latter dichotomy appears reminiscent of its well-known non-quantum real- 
space analogue, but the parallel is misleading because in quantum case the freedom 
of employing an additional coupling constant 7 enables us to circumvent the re- 
strictions. Indeed, once we select any dimension D, angular-momentum index m, 
winding number N and any "allowed" integer M, our spectral recipe (jTSj) may sim- 
ply be re-read as an explicit definition of the knot-supporting value of the coupling 
constant 

[M\ 2 ( D-2\ 2 

7 = bivJ -\ m+ —) ■ 

This implies that at non-vanishing 7s, the quantum knots do exist in any dimension. 

5 Discussion 

In the language of physics, our present construction and solution of a new and fairly 
unusual exactly solvable quantum model of bound states is based on the freedom 
of choosing the knot-shaped, complex contours of integration C. This trick is not 
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new [2J [3] and may be perceived as just a consequence of the admitted loss of the 
observability of the coordinates in PT-symmetric Quantum Mechanics [18J. 

From an experimentalist's point of view, the omission of the standard assumption 
that the coordinate "should be" an observable quantity is not entirely unacceptable 
since the current use of the concept of quasi-particles paved the way for similar 
constructions. Related Hamiltonians could be called, in certain sense, manifestly 
non-Hermitian. Still, they are currently finding applications in nuclear physics (where 
they are called quasi- Hermit ian [TO])- The loss of the reality of the coordinates is 
also quite common in field theory where the similar unusual Hamiltonians are being 
rather called CPT-symmetric [19] or crypto- Hermitian [20] . 

In a pragmatic phenomenological setting, the fairly unusual nature of the new 
structures of spectra seems promising. At the same time, the formalism itself is now 
considered fully consistent with the standard postulates of quantum theory. In the 
language of mathematics, the emergence of its innovative features may be understood 
as related to non-locality, i.e., to the replacement of the standard scalar product 



Although this leaves an overall mathematical consistency and physical theoretical 
framework of Quantum Theory virtually unchanged [18] , a new space is being open, 
inter alii, to the topology-based innovations. In principle, they might inspire new 
developments of some of the older successful applications of the formalism ranging 
from innovative supersymmetric constructions [22] to cosmology [23], occasionally 
even leaving the domain of quantum physics [24]. 

In section [2] we summarized briefly the key ingredients of quantum theory where 
the "correct" metric is assumed nontrivial, O^J. Let us now add a few comments 
which may have emerged during our subsequent transition to the quantum-knot 
models of section |4J Of course, our eq. ([2]) at N > 1 can still be treated as compatible 





by its generalized, nonlocal modifications [HI [21] 





13 



with the standard postulates of quantum theory in principle. We only have to repeat 
that the necessary proof of the latter compatibility statement is nontrivial. For each 
individual Hamiltonian (with the property H ^ W with respect to the specific 
Dirac's definition of the t— conjugation) the rigorous demonstration is indispensable 
that the spectra are real and that they are discrete and bounded from below. This 
demonstration represents, in fact, the main part of our present contribution. 

It is precisely the difficulty of the latter step which motivated our present start 
from the dynamically trivial version (j2J) of the BG model with vanishing A. In the 
nearest future we shall have to pay attention to the related operator Q(x,y), feeling 
inspired by the Mostafazadeh's [25J explicit formula 

M 

e(x,y) « £ y n (x)s n V n (y), M>1 (19) 

n=0 

where the normalized eigenstates ^ n {x) of W have to be constructed in ?^( ori 9 mai ) 
and where the real and positive constants s n > 1 are, in principle, arbitrary [26J. 
Of course, the letter-format of our present message does not allow us to get too 
far beyond the citation of the encouraging observation that this type of formula 
exhibited a quick convergence to the exact in the square-well model where a fairly 
good approximation has already been obtained at M ~ 10 [I]. 

6 Summary 

In spite of the absence of any confining force, our Schrodinger eq. (j2j) defined along 
topologically nontrivial integration paths has been shown to generate certain bound 
states ip(r) G L 2 (c^^j at a discrete set of the centrifugal coupling 7. We may 
emphasize that in such an exemplification of the more or less standard quantum 
theory 

• the complexified coordinates are loosing their immediate observability, 

• arbitrary complex potentials V(r) eC are allowed, provided only that the 
spectrum remains real, 
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• a redefinition of the inner product in the Hilbert space is required in order to 
return to the standard probabilistic framework of quantum theory, 

• a challenging general open problem arises concerning the role and tractability 
of the complex coordinate paths with a nontrivial topological structure. 

Of course, our present, exactly solvable N > 1 quantum- knot bound-state problem 
would probably become purely numerical after its immersion in virtually any external 
confining potential. In this setting, even the question of survival of the reality of the 
new bound-state spectra at A 7^ remains open. 
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Figure captions 

Figure 1. Sample of the curve with N = 1. 

Figure 2. Sample of the curve with N = 2. 

Figure 3. Sample of the curve with N = 3. 
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